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CN ■ I. INTRODUCTION 

In the context of Quantum Field Theory (QFT) a rich non-perturbative vacuum structure associated with the 
mixing of fermion and boson fields has been revealed [1-15] and the exact formulas for fermion and boson field 
^ ' oscillations [2,7,10-12] are now established. In particular, in a full QFT treatment particle mixing exhibits new 
features with respect to the usual formulae in Quantum Mechanics QM [16]. The phenomenological analysis for 
^vq meson mixing has shown that, while for most of mixed systems the non-perturbative structure of the vacuum produces 
^ ' negligibly small effects, for strongly mixed systems such as co — (j) or rj — rj ' non-perturbative corrections can be as 
00 large as 5 — 20% [14]. Moreover the non-perturbative field theory effects may contribute in a crucial way in other 
physical contexts. For example, as shown in [17], the neutrino mixing may contribute to the value of the cosmological 
constant exactly because of the non-perturbative effects. 

There are, however, several aspects which still need to be fully developed. For example, how to deal, in the presence 
[ of mixing, with those decay processes where neutrinos are generated. 
' Since the time of the introduction of the Pontecorvo mixing transformations [16], it is well known that the mixed 
(flavor) neutrinos are not mass eigenstates. This implies that flavor neutrinos are not representations of the Poincare 
group: one cannot think of them as asymptotic fields in the frame of the Lehmann-Symanzik-Zimmerman (LSZ) 
formalism [18]. The QFT analysis of the mixing phenomenon has indeed clarified [1] that flavor neutrino field operators 
Oh, do not have the mathematical characterization necessary in order to be defined as asymptotic field operators acting 
^ ' on the massive neutrino vacuum. The origin of this is related to the fact that the vacuum for the massive neutrinos 
turns out to be unitarily inequivalent to the vacuum for the mixed neutrino fields. 

Previous works were mainly focused on the determination of the oscillation probability through the analysis of the 
expectation value of the flavor charge operator 
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X 

Wpit)\QMt)) (1) 

Q being a function of the flavor annihilation/creation operators [19] and p,a — e,fi,T flavor indices. In fact, flavor 
states are produced in weak interaction processes and we are left to the question whether Eq.(l) can be consistently 
extended to include the neutrino production vertex and what is the explicit form of 

(n(t)|Q,jnw) (2) 

where |n(t)) represents the evolution at time t of the parent state for the neutrino. 

The computation of the matrix element (2) is not trivial matter since it simultaneously involves LSZ states and flavor 
states, which are not LSZ. This paper is devoted to the study of such matrix element and of several related topics. 
In Section II we review the formalism of neutrino mixing in QFT and discuss the proper deflnition of flavor charges 
and states; in Section III we clarify the issues of the possibility of different mass parameters in field mixing, which 
has recently attracted some attention [5,7,9,10,20]. In Section IV we perform a careful analysis of (2). Conclusion are 
drawn in Section V. For the reader convenience, we present in Appendix A the proof of orthogonality of the flavor 
states at different times. The Appendix B contains the explicit form of some equations derived in the text. 
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II. FLAVOR FIELDS AND STATES 



Let us start by introducing the general frame for our discussion, which is also useful to set up our notations. 
For a detailed review see [13]. For simplicity we consider only two Dirac neutrino fields. The Pontecorvo mixing 
transformations are [16] 



Ve {x) = Vi [x) COS 9 + V2 {x) sin 6 
>^ti{x) = —i'i{x) sin 9 + i^2{x) cos 9 , 



(3) 



where 9 is the mixing angle and Pi and 1/2 are massive non-interacting, free fields, anticommuting with each other at 
any space-time point. The fields vi and U2 have non-zero masses mi ^ TO2 and are explicitly given by 



1,2. 



(4) 



with = e~*'^'=''*uj^ j(0), vl^ ^{t) = e*'^'''**i'k^j(0) and ujk,i = + mf. Here and in the following we use t = xq, 

when no misunderstanding arises. The vacuum for the ai and (3i operators is denoted by |0)i,2: CKki|0)i2 = 
/JJ^ J0)i2 = 0. The anticommutation relations are the usual ones (see Ref. [1]). The orthonormality and completeness 
relations are: 



(5) 



The fields and are completely determined through Eq.(3), which can be rewritten in the form (we use 
(a,i) = (e,l),(M,2)): 



(6) 



with Geit) the generator of the mixing transformations Eq.(3) 



Ge{t) = exp 



9 j (fx (^i'l{x)i'2{x) - v'l{x)ui{xfj 



(7) 



Eq.(6) provides an expansion of the flavor fields Vg and V/^ in the same basis of Ui and i/2. The flavor annihilation 
operators are then identified with 



(8) 



The flavor vacuum is defined as \0{t))e^fj, = G'^^(t)|0)i.2 and turns out to be orthogonal to the vacimm for the mass 
eigenstates |0)i_2 in the infinite volume limit. Note the time dependence of \0{t))e,fi'. it turns out that flavor vacua 
taken at different times are orthogonal in the infinite volume limit (see Appendix A). In the following for simplicity 
we will use the notation \0)e,^c = |0(0))e,;^ to denote the flavor vacuum state at the reference time t = 0. 
The explicit expression of the flavor annihilation/creation operators for k = (0, 0, |k|) is: 



ce 







seVk(t) 



-seU^{t) ce seV^{t) 

-50 14(0 ce S0U*{t) 
-seVk{t) -seUk{t) ce 



where ce = cos 9, se = sin 9 and 



r 



(9) 



(10) 
(11) 



with = (-1)''. We have: 
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V 2uJk,i J \ 2cjfe,2 / V {^k,i + mi){u>k,2 + ■m2) J 

,y>_( '^k,i+mi \^/ u!k,2 + m2 \^/ |k| |k| \ 

^ \ 2u)k,\ ) \ 2ujk,2 ) \(wfe,2+m2) (wfc,i+mi)y' 

|f/k|' + = 1. (15) 

As discussed in Ref. [3], in the two flavor mixing case the group structure associated with mixing transformations 
is SU{2) and one can define the following charges in the mass basis: 



i{t) = \ j S^ul{x)TjVm{x), j = 1,2,3, (16) 



where = {1^1,1^2) and Tj = (Jj/2 with being the Pauli matrices. The U{1) Noether charges associated with vi 
and V2 can be then expressed as 

Ql = + Qm,3 ; Q2 = - Qm,3- (17) 

with Q total (conserved) charge. As usual, we need to normal order such charge operators: 

:Qi-= Jd'^: ,yl{x) u,{x) := J d'i^ ' ^Xil^-Ki) ' * = 1' 2. (18) 

where the : .. : denotes normal ordering with respect to the vacuum 1 0)1,2. 
It is then clear that the neutrino states with definite masses defined as 

kk,i)="k>)i,2, i=l,2, (19) 

are eigenstates of the above conserved charges, which can be identified with the lepton charges in the absence of 
mixing. 

The situation changes when wc turn to the flavor basis [3]: 

QfAt) = \j <f^v\{x)TjVf{x), f = e,M, i = l,2 (20) 

where vj = [u^^Uij). Observe that the diagonal SU{2) generator Qf,i{t) is time-dependent in the flavor basis. 
Thus the flavor charges defined as 

QuAt) = \Q + Q}At) ; Q-.W = ^0-Q/,3(0, (21) 

arc now timc-dopcndcnt and are the lepton charges in presence of mixing [3]. Indeed their expectation values in the 
flavor state (see Eq.(26) below) give the oscillation formulas [2,7,10-12]. 

Particular attention has to be paid now to the normal ordering issue. We deflne the normal ordered charges 
:: Qv„{t) :: with respect to the vacuum \0)e,ij. as 

::g.^(i)::^ J d'^ :: ut{x) Mx) :: = ^ J d^k (^a'^l^{t)al^^^^^^ a = e,n;, (22) 

where the new symbol :: ... :: for the normal ordering was introduced to remember that it refers to the flavor vacuum. 
The deflnition for any operator A, is the following 

:: A::= A - e,^(0m0)e,^ (23) 

Note that: 
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" Q^At) :: = Gg\t) : Qj : Geit) , with {a,j) = (e, 1), (/x, 2), (24) 

and 

:: QuAt) " + " QuAt) ■■ = ■■ Qi ■■ + ■■ Q2 ■■ = ■■ Q ■■ ■ (25) 

We define the flavor states as eigenstates of the flavor charges Q^^ at a reference time t = 0: 

= <lA^)\m)e,^, CT = e,ll (26) 

and similar ones for antiparticles. We have 

QuAO) l^'k.e) = kk,e) ; Q..(0) Wl,^) = Wl,^) (27) 

and :: Q.JO) :: 1^.'^,^) =- Q.,(0) :: = 0- Moreover 

"QuAQ) ■■■■ |0)e,p = 0. (28) 
These results are far from being trivial since the usual Pontecorvo states [16] : 

= cose + sine |i/'k,2) (29) 

kk,^)p = - sine + cos^ \ul^2) , (30) 

are not eigenstates of the flavor charges, as can be easily checked. 

It is instructive to consider the expectation values of the flavor charges onto the Pontecorvo states, in order to 
better appreciate how much the lepton charge is violated in the usual quantum mechanical states. We find: 

pK,el "<3-e(0) " WlJp = cos^0 + sm'^e + 2\U),\sm'^ecos'^e + J2 [d% (31) 

and 

i,2(0|::g.,(0)::|0)i,2 = 5^ / d^K (32) 

r •' 

Eqs.(31) and (32) clearly are both infinite. 

One may think that the problem with infinity is due to the normal ordering with respect to the flavor vacuum and 
consider the cxp{!ctation values of : Qv^{t) :, i.e. the normal ordered flavor charges with respect to the mass vacuum 
|0)i,2. One has then 

pK,eh <3^e(0):|z^k,e)p = cos* 6' + sin* 6' + 2|C/k|sin2 6»cos^ 61 <1, V6i ^ 0, mi ^ ms, k ^ 0, (33) 

i,2(0| :Q.,(0) : |0)i,2 =0, (34) 



and 



i,2(0|(: g.,(0) :)2|0)i,2 =4sin2ecos2e J d^k\^,f, (35) 



pK,el(: Q-^A^) ■.fWl,e)p = cos'^0 + sin^0 + sin^0cos'^e 



2|«7k| + |C/k|'+4 / d^klVi 



(36) 



which are both inflnitc, thus making the corresponding quantum fluctuations divergent. 

Hence, we conclude that the correct flavor state and normal ordered operators are those deflned in Eqs.(26) and 
(23) respectively. 
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III. MASS PARAMETERS AND FIELD MIXING 



In Eq.(6) and v^j^ ^ are the spinor wavefunctions of the massive neutrinos Vi, i = 1,2. As already observed in 
the previous section, Eq.(6) provides an expansion of the flavor fields i^o-, cr = e, /i, in the same basis of i = 1,2. 
However, it was noticed in Ref. [5] that expanding the flavor fields in the same basis as the (free) fields with definite 
masses is actually a special choice, and that in principle a more general possibility exists. Indeed, in the expansion 
(6) one could use cigcnfunctions with arbitrary masses fi^- In other words, the transformation (8) can be generalized 
by writing the flavor fields as [5] 

u,{x) = E [<AtWu,<^. it) + ^-k,.(i)^:l,.„ it)] e*•^ (37) 

where u„ and Va- are the helicity eigenfunctions with mass /i^ (the use of such helicity eigenfunctions as a basis 
simplifies calculation with respect to the choice of Ref. [1]). In Eq.(37) the generalized flavor operators are denoted 
by a tilde in order to distinguish them from the ones defined in Eq.(8). The expansion Eq.(37) is more general than 
the one in Eq.(6) since the latter corresponds to the particular choice ji^ = m\, = m^- 

Since the issue of the arbitrary mass parameters in the field mixing formalism has attracted some attention 
[5,7,9,10,20], it is worth clarifying some basic facts about the choice of mass parameters within QFT in general, 
independently from the occurrence of the field mixing phenomenon. 

We will refer to fermion fields since in this paper we are interested in neutrinos, but the conclusions can be also 
extended to boson fields. 

First of all, it is worth noting that the mass parametrization problem can be revealed also in the free field case. 
Indeed, one may still consider the change of mass parametrization m — > /U, which correspond to choosing uj^, wj^ as 
free field amplitudes with the arbitrary mass parameter fi. 

Consider the set of free (fermion) field operators composed, for simplicity, by only two elements, i.e. assume our 
operators are 

"t'* ) , with i = 1,2. (38) 

Let the non-zero masses be mi and m2, with mi ^ m2. The wave functions Wj and Vi (we omit the index k whenever 
no confusion arises) satisfy the free Dirac equations 

mi)ui = , {i}t- mi)vi = , (39) 

respectively. Let |0)i^2 be the vacuum state annihilated by ■ and p'iy^ ^. 

Since in QFT there exist infinitely many unitarily inequivalent representations of the canonical (anti-)commutation 
relations [18,21], one could consider another fermion set of operators, say 

related to Eq.(38) by a Bogoliubov transformation (see below). The freedom of choosing another set of operators is, 
for example, typically exploited in QFT at finite temperature, or more generally when one introduces the irreducible 
set of "bare" field operators in terms of which the Lagrangian of the theory is written. In such a case the set of bare 
fields is not necessarily composed by the same number of elements as the one of the set of physical (asymptotic) fields 
satisfying free field equations and in terms of which observables are expressed. In general, indeed, bound states of 
bare fields may also belong to the set of physical fields. The mapping between the bare fields and the asymptotic 
fields is called the Haag expansion [18,21]. 

Suppose the wave functions of the field operators in Eq.(40) also satisfy the free Dirac equations 

+ Hi)ui = , injui = , (41) 

respectively. The mass parameter /ij in Eqs.(41) represents now the mass of the corresponding arbitrarily chosen 
fields in Eq.(40) and therefore it represents an arbitrary parameter. Our two set of operators are related by the 
transformation 
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with 



(43) 



where = (xi — Xi)/2 and cotxi = |k|//Zi, cotXi = |k|/m,. Notice that for /ii = mi, /i2 = ?7i2 one has I^{t) = 1, as 
it must be for the identity transformation. 

The explicit matrix form of Eq.(42), written for both i values, is: 





\ 





Ptit) 

















( ^ 




ak,2 











(44) 



where 



COS^j 5r 



-i{uik,i+Ulk,i) 



with i = 1,2 and LOk,i = -^/k^ + /x?. The vacuum state annihilated by the {al^^-, PL^^i) operators is 

|o(i))i,2 = /;'(t)|o)i,2 . 



(45) 
(46) 



(47) 



We observe that Eq.(42) is indeed nothing but the Bogoliubov transformation which relates the field operators {a^ ^, 

f^-ki) and (a{^j, /3!Lki)> of masses rrii and fii, respectively. In the infinite volmiie limit, the Hilbert spaces where 
the operators a^. and aj^ are respectively defined turn out to be unitarily incquivalent spaces. Moreover, the 
transformation parameter acts as a label specifying Hilbert spaces unitarily incquivalent among themselves for 
each (different) value of the /Xj mass parameter. 

We note that the vacuum |0(t))i.2 is not annihilated by ^ and ^ and it is not eigenstate of the the number oper- 
ators = 5^,. / d'^k a^^al- ^ and Np. = J (P'k /3^^j/3^ j. Similarly |0)i.2 is not annihilated by ,■ and f3^ ^ and it 

is not eigenstate of the the number operators Na^it) = J2r I ^'^^ ii^)^k ii^) ^'^d Np^{t) = J2r I (^kii^)(^L ii'^)- 
One obtains 



1,2 



2tk 



(0(t)|iV„jO(t))i,2 = i,2(0(t)|iV^jO(t))i,2 = sin^e 



and 



2ck 



i,2(0|iVa,(0|0)i,2 = i,2(0|7V^,(t)|0)i,2 = sin^e 



(48) 



(49) 



In other words, the number operator, say Na-, is not an invariant quantity imder the Bogoliubov transformation 
Eq.(42): it gets a dependence on the mass parameters. This is, however, not surprising since it is known that the 
Bogoliubov transformation Eq.(42) introduces a new set of canonical operators and a new (i.e. in the infinite volume 
limit unitarily, and therefore physically, incquivalent) Hilbert spae(\ Stated differently, through the Bogoliubov 
transformation a new set of asymptotic fields (a new set of quasipartic:les) is introduced, i.e. there are infinitely many 
sets of asymptotic fields, each set being associated to its specific representation. The choice of which one is the set to 
be used is then dictated by the physical conditions which are actually realized. For example, the mass values which 
have to be singled out in the renormalization procedure must be the observed physical masses. 

Since the tilde quantities correspond to some new quasi-particlc objects and the tilde nmnber operator describes a 
different type of particles, then, the number operator average shall not be expected to remain the same under such 
transformations. Indeed, defined the state |5j^ J as |5j^ ^(0)) 



Q:^^-(0)|0), we have: 



(51c,i(0) I N^M I 51,,,(0)) = \{al,M,^l.limf = ||p,^|2e*(-M-M)* + |^k|2g,(^,,.+..,.)t|2^ 



i = l,2. 



(50) 



which shows that the expectation value of the time dependent number operator is not preserved by the transformation 
(42) applied to both states and operators. Nevertheless, in the cases of free fields, the charge operator is still conserved 
in transformation like Eq.(42): 
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(51) 



moreover the expectation value of the charge at time t on the state at time t = is also free from mass parameters 

{alM I Qi{t) I al^.m = {alM \ Q,(i) | ai;,,(0)), (52) 

since we have 

WuM^ilim' + mlM^Zmi' = i{<w,<.(o)}p + \{PilM,<^m\'- (53) 

Similar results are obtained in the two flavor particle mixing case. Note, however, that in the case of three flavor 
mixing and in presence of CP violation, the charge operator and the flavor states are dependent on the arbitrary mass 
parameters and the quantity which is invariant under the transformation like the Eq.(42) is 

«.J l"k,.J - e,^,r(0| :: Q.„(i) :: |0)e,^,., (54) 

with |aky_^) = (0)|0)e,/i,T and cr = e,/U, r as observed in Ref. [10] 

Having clarified that the possibility of difii'erent mass parameters is intrinsic to the very same structure of QFT and 
is independent of the occurrence or not of the field mixing, we may aflirm that the mass parameters must be chosen 
not arbitrarily, but they must be justified on the ground of physical reasons [4,9,10]. 

In particular, in the mixing problem, the choice /ig = mi, /i^ = m2, /ir = rn^, is motivated by the fact that mi, 
1712 and ma are the masses of the mass eigenfields and therefore such a choice is the only one physically relevant. 

In our computations, instead of using the number operator, we use the charge operator Q which in the mixing 
phenomena describes the relative population densities of flavor particles in the beam and it is related with the 
oscillating observables that are; the lepton charge, in the case of neutrino mixing, the strange charge in meson 
systems like — and B'^ — B'^, the charmed charge in the systems — D^, and so on. We note that, as shown in 
refs. [12,14], the momentum operator for the mixing of neutral fields plays an analogous role to the one of the charge 
operator for charged fields. 



IV. EXPECTATION VALUE OF NEUTRINO LEPTON CHARGE IN THE ONE PION STATE IN PION 

DECAY 

In this Section we study the structure of the flavor charge expectation values in the case the production process of 
neutrinos is taken into account. This is done by using the flavor Hilbcrt space discussed above. In Ref. [22] a similar 
calculation was performed by using the mass Hilbert space, so neglecting flavor vacuum effects. 

The flnal aim of the authors of Ref. [22] was to derive an oscillation formula in space, which is relevant for current 
experiments. On the other hand, a general oscillation formula with space-time dependence has been obtained in Ref. 
[11] in terms of expectation values of the flavor currents on the flavor neutrino states, exhibiting the corrections due 
to the flavor vacuum. In the following, we show in an explicit way how calculations can be performed with interacting 
fields on the flavor Hilbert space. 

We consider the case where the flavor neutrinos are produced through the pion decay 7r+ — > yf^ +v^. We use the 
phenomenological approach to the pion decay [19] without referring to the quark structure of the pion. As initial 7r+ 
state at time we use ]n((k),x, a;^)) = |7r((k),x, a;^)) x \Q{x^j))„ x [0)^ with (k) the average k vector (see below). 



The neutrino vacuum is the flavor vacuum: |0(a:j))^ = Gq ^(xj)|0)i,2 (note the change of notation for |0(a;°)) 
with respect to the previous sections). We calculate the expectation values of neutrino flavor charge :: Qvu,{x^) " with 
respect to one pion state |H((k),x, a:°)) where < .x" 

(:: Q.^ix') ::) = (n((k), x, x")] :: Q,^{x°) :: |n((k), x, a:")) 

= (n((k),x,x?)|5-^(xO,x?) :: g.^(x°) :: S {x° , x°)\U{{k) , x°)) (55) 

with 

|n((k),x,a;0)) = 5(a;0,ar?)|n((k),x,a;?)) (56) 
where, in the interaction representation, the normal ordered neutrino flavor charge :: Qi,^{x^) :: is 
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and S{x°,x°) is 



;(ym)- 



The muon neutrino field and the muon field are expanded, in the interaction representation, as 



„ J (ZTt) 2 1" 



.0\/3rt ^^0^„-iq•x 



(58) 



(59) 



(60) 



with Q:q_^^(a;°) flavor annihilation operator for neutrino. In Eq.(59), u^^2{^^) is the spinor wave function for the 
massive neutrino field 1^2(2;). 

The hamiltonian of weak interaction is 



i?int(a;) = -i Qt, Vij,{x) 7^(1 + 7'"') /i(.T) dxTr{x) + h.c, 

where 7r(a;), /u(x), f^(x), are the fields of pion, muon and fiavor (muon) neutrino, respectively. 
In the lowest order (i.e. the second order) of the weak interaction, we have 



(61) 



(:: Q.^(xO) ::) = (n((k),x,x?)| 



/ rfSFi„t(^) "Q.^(a;°) :: / (fyH,, 



|n((k),x,x?)) 



and, in the case of positively charged pion 7r+ , we have explicitly 



(::Q.^(x°)::)= / / ^^^^(^((k), x, a;°)| 

Jx» Jx° 



(62) 



(63) 



X [igldxn\z) ^(0|77(-2)7^(1 +75) .(0(x?)|zv^(^) :: g.^(x°) :: I7^(y)|0(x?)). 
X 7<^(1 + 75)M(y)|0)^ «57raa7r(y)] |7r((k), x, x?)). 
Expressing the charge operator as in Eq.(57), a straightforward calculation gives 

(:: Q,^{x^) ::) = / / rf*^A(^((k), x, ar?)|aA7rt(^)a,7r(y)|7r((k), x, x?)) 

^ ^ d^p rf3q ^3q/ ^3q,, ^^.(^^^^..^..(p+q,).^ 



r,t,u,v 



X e-'(-^+--^)-°e'(--+-.'.^)^''z;^,^7^1 + 75)<,2"q',27"(l + 75)'^';,,^ 
+ similar terms, 



(64) 



being uip = Jp'^ + mi and ujq^2 = V^?-H 



The explicit expression of Eq.(64) is given by Eq.(Bl) in Appendix B. 

Evaluation of the Bosonic term 

We consider first the bosonic part: {7r{{k),x,x^)\d\'jT^{z)da'JT{y)\Tr{{k),x,x^)) . The pion state is defined as: 

^^(k, (k)) 6-*-^+'"'=^? 



|7r((k),x,ar?)) = J d?k 



(65) 
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with 

\2 



and Ofc = ^J\s?' + m%. The pion field in the interaction picture is: 

<x) = J (07^ {^^ + bl e-C' -^^^*)) , (67) 

hence 

(7r((k),x, x°j)\d^nHz)dMy)Wm,^, ^?)) = 



/: 



4(27r)3f2fcf2fc 
thus Eq.(64) may be expressed as 

X g-i(^-P-^)-^gi(^' -P-'i')-y gi(^l'-i<'p-uiq,2){z° -x^j) ^-i{n^, -U!p-u>^, ^2){y° -^°) y'^ ^7'^(1 + 75)Wq 2Wq' 2 

X 7'^(1 + 75)^;^,^ .(0(a;?)|<,.^(^°)[<t ,^(a;0)^.,,^^^(^0) _ ^.t 

+ similar terms, (69) 
(see Eq.(B2) in Appendix B) and then 

d^k d^k' d^p d\ /^,^^A_(^.' /w\w^(k-k')-Xe-i(n.-a,,)x? 



X .(0(x°)|aLp,.,(^°)Kl.,(^°)<,.,(^°) - /?qk(^°)/3q,.,(^°)]«k?-p,.,(j/°)|0(x?)). 

+ similar terms, (70) 

(see Eq.(B3) in Appendix B). 
Evaluation of terms like 

.(0(x?) |aL-p,., (^°) - /3q!., (^°)] «k^-p,., (2/°) |0(x?)). 

We have 

. (0(x?)|aL_p,.^(z")K|.^(x"X,,^(x°) - /3;;|.,(x°)/3:;,.^(x")]a;:Lp,.^(2/°)|0(x°)). = 

= i,2(0|G(a.?)4-p,.,(^°)G-^(a;?)G(a;?)[<t^^(a;0)<,,^(a.O) 
-/3^!.,(a=°)/3^,.,(a;°)]G-i(a.?)G(a;?)<t_p,.,(y°)G-i(a;?)|0)i,2 

= iMO\&L-p,uyj,z°)[all^ix°,x')a^^^^^^^ (71) 
Where the flavor annihilation operators are a^^,^{t,6) = Gg^{t) a^^^ G0{t): 
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<..(i,^) = cos^a^,i + sin^ ^ [w;|i(t)u^,2(i)<,2 + <!iW^^lq,2(i) ^l^q,2 

S 

a;,,^(t,0)=cos0a^,2 - sine ^ + ^sW^lq.iW /^-q,i 

/3!:q,.e(*.^) = cose/5-q,i + Sine + <!2(iKq,i(*) "qV 

/3-q,.,(i,^)=COSe/31q,2 - sine J2[^-c,,lit>-^,2{*) + <!lW^lq,2W <!l 

By using d^,,„(t',t) = Ge(t') "q.^J*,^) Gg\t'), and G^'W = ^-^(0, we have 

+ sine + K\^{t)v\2{t)P'\,.,{t',-e) 

S 

d^,,^(i',i) = cose<^,^(t',-e) 

s 

Cq,.e(*',*) = cose/3:q,,^(t',-e) 

/31q,.,(i',i) = cose/31q^,^(t',-e) 

- sine ^[^i;,i(i)^:q,2W/?lq,.e(*'.-^) + W^lq,2 (t' , "^^ 



(72) 



(73) 



By using Eqs.(72) in Eqs.(73) the last equality in Eqs.(71) can be finally expressed in terms of the massive neutrino 
operators a^^^, a^'t, /Jllq^j, P-^^i which act on the massive neutrino vacuum |0)i,2- Thus our computation allows 
the expression of the matrix element {Il{t)\Qi,^\U{t)) as a function of states that have the proper mathematical 
characterization of the LSZ formalism. 

We observe that the formula (69) is consistent with Eq.(14) of Ref. [22]: indeed, when the pion state is represented 
by a plane wave (as done in Ref. [22]), the two equations acquire the same form, except for the neutrino charge 
expectation value, which in our case includes the flavor vacuum contributions. If we neglect such flavor vacuum effect, 
we precisely recover Eq.(14) of Ref. [22]. 

Furthermore, when we consider the effect of flnite lifetime of the pion and we take the neutrino oscillation time 
— x'j much bigger than l/F^, then a^^ {y'^)\0(x'j))i, ~ W\i f^i^'j)) wc obtain a form which is similar (once 
integrated in time) to the one considered in Ref. [11] giving the space dependent oscillation formula, except for the 
fact that the pion state rather than the neutrino state is represented by a wave packet. The flnal expression ensuing 
the formula (69) is too lengthy to be discussed at this stage. Further analysis will be done elsewhere. 



V. CONCLUSION 



In this paper wc have studied neutrino mixing and oscillations in quantum field theory and we discussed the 
determination of the oscillation probability including the neutrino production vertex. A crucial point in our analysis 
is the disclosure of the fact that in order to describe the neutrino oscillations we have to use the flavor states deflned 
as |i^^ct) ~ '^kv (0)|0(0))i/, with a = cJ-I.t and the flavor charge operators :: Q,^^ :: normal ordered with respect to 
the flavor vacuum. Indeed, we have shown that the usual Pontecorvo states are not eigenstates of the flavor charges 
: : and :: ::. 

Wc showed that the possibility of different mass parameters is intrinsic to the very same structure of QFT and is 
independent of the occurrence or not of the fleld mixing; hence, as noted in [4,9,10], the mass parameters must be 
chosen not arbitrarily, but on the ground of physical reasons: /Xg = rni, /i^ = m2, /Xt = n^s- 

Moreover, we have computed the neutrino Icpton charge in decay processes where neutrinos arc generated, proving 
that the corresponding lepton charge expectation value can be uniquely expressed in terms of LSZ states and in 
particular of the massive neutrino annihilation/creation operators acting on the massive neutrino vacuum. 
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APPENDIX A: ORTHOGONALITY OF FLAVOR STATES AT DIFFERENT TIMES 

The product of two vacuum states at different times t^t' (we put for simplicity t' = 0) is 

.(0|0(t)), =l[Cl{t) = e^S^'"^''^ (Al) 

k 

with 

Ck{t) = (l-sin^^ 2 sin^^ cos'' 6 \Vk\^ e-iio^k,2+'^k,i)t 

In the infinite volume limit we obtain (note that |Cfc(/:)| < 1 for any value of k, t, and of the parameters 0, mi, m2 ): 

2V 



lim i/(0|0(t))^ = lim exp 

V — ^oo V — ^oo 



(27r)3 



d^k {In \Ck{t)\ + iak{t)) 



(A3) 



with \Ck{t)\' = i?e[Cfc(t)]2 + /m[Cfc(t)]2 and ak{t) = tan"! {Im[Ck{t)]/Re[Ck{t)]). 
Thus we have orthogonality of the vacua at different times^ . 

One can easily check that flavor neutrino states are also orthogonal at different times. Consider the electron neutrino 
state at time t with momentum k: 

l'^k,eW)=<..W|0(0).- (A4) 

The flavor vacuum is explicitly given by 

|0(t)). = nGp,Ui)|0)i,2, (A5) 
p 

where, to be precise, the mass vacuum is to be understood as |0)i,2 = |0)i 2 |0)i 2 |0)i 2--" 
Then, we have 

K,emKM = AQ\^i.soxlst)m). (A6) 

= (ni.2(0|Gp,.(0)j aL,..(0)<.,(<) (jl 0^,^010)1,2^ , (A7) 
and, since for p ^ q the mixing generators commute among themselves and with aj^ for k 7^ p, q, it is 

WL,emK,em ^ >'(0''K.e(0)<.eW|0''W>- n i,2(0|Gp,,(0)G-;,(t)|0)i,2 (A8) 

= AO''\aluSO)allSm''{t)). .(0|0(t)).. (A9) 

By using Eq.(A3), we obtain the orthogonality of flavor states at different times in the infinite volume limit, provided 
i^(0''|ajj. ^ {Q)oikiy {t)\0^{t))i' is finite or zero, as it is indeed. 



^Note that it may be that for some values of the continuous index k, the Ck{t) arc periodic functions of time. However, the 
set of values of k for which this happens is a zero- measure set in the above integration Eq.(A3). 
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APPENDIX B: EXPLICIT EXPRESSIONS OF EQS.(64), (69), (70) 



The complete Eq.(64) is: 

(::Q.^(a;0)::) = -|5.P / / rfSAW(k), x, x?)|aA7rt(^)a,7r(y)|7r((k), x, ar?)) 

r.t.u.v ^ ^^""-^ 



X .(0(a;?)|<,,^(^0)[<t_,^(:r°)<„,,^(x°) - /3;t,.^(a:°)/3^.,.^(x°)]a:;t,^(yO)|0(a;?)). 
+ e'(P-'i)-^e-'(P-'i')-ye-'('^''-'^<''='^°e'('"''-'^<''.^'f°v^_^7^(l + 75X,2^^q',27''(l + 75)i^p,^ 
X .(0(:c?)|/3^t^^(^0)[a;t,_,^(xO)a:;,, ,^(a;°) - /^^^^.^(x")/?;;.,.^^")]/?:;, ,JyO)|0(x?)). 
+ e'(P+'i>'=e-'(P-'»'>ye-*('^''+'^''-^)^°e'('"''-'^<''.^'^°v;;^^7^(l + 75X,2'^q',27''(l + 75)fp,^ 

+ e'(P-'»)-^e-'(P+'»'>ye-*('^''-'^<''=)^°e'('^''+'^<''.^'^%;,,^7^(l + 75)<,2«q',27"(l + 75)^'p,^ 

X .(0(x?)|/3*t^^(^0)[a^L^^^(a;")<.,.^(x") - /3;I,.^(xO)/3;i,, .^(x")]<t^(y°)|0(x?)).) (Bl) 

The complete Eq.(69) is: 

X ^e~*'''~P~''^''^e*''''~P~'''^'^e'^^''~'^''~'^'''^^^^°~^°^e~''^'^'~'^''~'^«''2)(y°-3J?)^r ^^-^(i -|- 75)^^ 2^^' 1 
X 7-(l +75)^;^_^ .(0(a;?)|<,,^(^°)[<t,,^(xO)<„,,^(a.O) - ,^(yO)|0(a;?)). 

^ g-»(k-p+q)-Zgi(k'-p+q')-yg»(f^fe-Wp+w<,,2)(z°-x?)g-i(n;,/-a>p+ai,/_2)(j/''-a;?)^-,r ^^^(X ^ 75)^^q 2''^q' 2 

X 7'^(1 +75)^;^,^ .(0(a;?)|/34t^^(^0)[^.t ^^^(^0^^.,, ^^(^0^ _ /3;t _^^(a;0)/3;;,, ,^(x°)]/3:;, ,^(y°)|0(a;?)). 

^ g-j(k-p-q)-Zgi(k'-p+qO-ygi(nfc-a;p-a;,,2)(z°-x?)g-i(nj,/-a;p+a;,/_2)(!/''-a:?),yr ^ 75)Wq 2'*^q' 2 

X 7'^(1 +75)^;^,^ .(0(a;?)|<,.^(z")[<t,,^(a:°)<„,,^(x°) - /3;t,.^(:r°)/3:;„,,^(x°)]/3;;,,,^(y°)|0(a;?)). 

^ g-i(k-p+q)-Zgi(k'-p-q')-ygi(nfc-a;p+w,,2)(z''-x?)g-i(afe/-Wp-a;,/_2)(!/''-a:?)^i- ^A^j^ _|_ -yg)^* ^u"/ 2 

X 7^(1 + 75)^p,^ .(0(x°)|/3*V,(^°)Kts.,(^°)<^.,(^°) - /3;t,,^(x°)/3^^^ 

(B2) 

The complete Eq.(70) is: 
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Jx° r,t,u,v'' \ ) k k' 

^ e^(^'=--''--''--=)(^"--?)e-'("'='--'•+-.-'=^=)(''°--?)^;;,^fc;,7^(l +75)4-p,2^^ fc;7"(l + 75)^^^,^ 

(B3) 
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